% tic;
% Objective
% f=[-3,-5,-3];
f=[-12;-9;0;0;0;0;];


% Inequality Constraints
% A=[];

% Equality Constraints
Aeq=[
    1 0 1 0 0 0
    0 1 0 1 0 0
    1 1 0 0 1 0
    4 2 0 0 0 1];
beq=[1000;1500;1750;4800];

% lower bound
lb=zeros(6,1);


% solver
% tic;
% x=linprog(f,A,b,Aeq,beq,lb,ub)
[x,fval,exitflag]=linprog(-f,A,b,Aeq,beq,lb)

% [x,fval]=linprog(f,A,b,Aeq,beq,lb,ub)
% [x,fval,exitflag,output, lambda]=linprog(f,A,b,Aeq,beq,lb,ub)
% toc;
